Lambda production in 40 AGeV/c Pb-Au collisions 
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Abstract. During the 1999 lead run, CERES has measured hadron and 
electron-pair production at 40 A GeV / c beam momentum with the spectrometer 
upgraded by the addition of a radial TPC. Here the analysis of A and A will be 
presented. 



1. Introduction 



CERES/NA45 is an experiment at the CERN SPS dedicated to measure low- mass 
e'^e" pairs near midrapidity in ultrarelativistic nuclear collisions [1]. The CERES 
spectrometer (Fig. 1) covers the full azimuthal acceptance in the polar angle region 
between 8° < 6* < 14.6°. The high electron identification capability allows to separate 
the rare leptonic signals from the large hadronic background; it is provided by two 
ring- imaging Cerenkov counters (RICH). Precise tracking of charged particles and 
vertex reconstruction are provided by two silicon drift detectors (SDD1,2) located 
behind a segmented Au target. This detector system is also used as a multiplicity 
trigger to get information on the centrality of the collision. 
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In order to improve the momentum resolution the CERES experiment has been 
upgraded [2, 3] during 1998 by the addition of a new magnet system and a Time 
Projection Chamber (TPC) with a radial electric drift field. The TPC is operated 
inside an inhomogeneous magnetic field with a maximal radial component of 0.5 T, and 
provides the measurement of up to twenty space points for each charged particle track. 
Besides a precise determination of the momentum the TPC also provides additional 
electron identification via dE/dx. The magnet system between the two RICHes was 
not operated in the upgraded configuration. 



TPC coils 
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Figure 1. The upgraded CERES spectrometer at the CERN SPS. In addition 
to the outlines of the various detectors are shown in the TPC field lines for the 
electric and magnetic field. 

The addition of the TPC substantially improves the hadron capability of the 
CERES spectrometer allowing a systematic investigation of hadronic observables 
around midrapity. In the fall of 1999 the upgraded CERES spectrometer was operated 
for the first time and 8 • 10^ semicentral Pb+Au collisions at 40 AGeV/c beam 
momentum were recorded. Because of a not yet completely functional read-out system 
this data set is limited in terms of statistics and momentum resolution. In the fall 
of 2000, the CERES spectrometer was fully operational with a very good overall 
performance, and 33 • 10® central 158 AGeV/c Pb+Au events were recorded. In this 
paper we will focus on the A analysis of the 40 AGeV/c data set. 

2. Analysis of the A hyperon signal 

In this A analysis we have used only the information from the SDD's and the TPC. 
The A hyperon and its antiparticle were identified by reconstructing their decays into 
final states containing only charged particles: A — > p 7r~ and A — > p tt"*". The 
geometrical acceptance for the A measurement is shown on the left side in Fig. 2. 
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We have also analyzed negatively charged particles (h~) and proton-like positive net 
charges, shortly denoted as (-f) — (— ) [4]. 
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Figure 2. Simulated geometrical acceptance of the TPC for A hyperons decaying 
into a proton and a pion. Left: no cuts, Middle: all cuts applied (see text). The 
box indicates the region selected by pair cuts. Right: Reconstructed decay point 
Vrec vs true decay point V;, in beam direction from a GEANT simulation. 



The momentum of charged particles is determined from their curvature in the 
TPC. Since we are not applying any particle identification method in the extraction of 
the strange particle signals we run through the combinations of all positive tracks with 
all negative tracks measured in the TPC assuming a decay hypothesis and calculate 
the invariant mass spectrum. In the same manner, the combinatorial background is 
evaluated using the mixed event method. 

To get a significant A signal we had to reduce the large combinatorial background 
using different types of cuts discussed in the following. 

Discrimination from target tracks: Taking into account the proper decay 
length ctq = 7.89 cm of the A and the 7-factor of about 5 for the TPC acceptance, 
the decay length in the laboratory is j3cTiab = P'ycTo ~ 40 cm. This means that more 
than 70 % of all A produced in the target decay after the SDD's (the distance between 
target and SDD2 is about 13 cm). To enhance the signal-to-background ratio, an 
optional cut requiring the absence of matching tracks in the SDD is used. 

Reconstruction of secondary vertex: In a second step, the decay vertex 
Vrec in z-direction (beam direction) is reconstructed for all A candidates by back- 
extrapolation of the TPC tracks and computing the point of their closest approach. 
This procedure was tested using a simulation. The result is shown in Fig. 2, right. 
This reconstructed secondary vertex is used in the analysis for an additional (optional) 
cut Vrec > 50 cm. At this threshold, the resolution in Vrec is 12 cm. 

Kinematical cuts: The extraction of the strange-particle signal is performed 
using geometrical and kinematical constraints. A fiducial cut of 130 < 9 < 240 mrad 
was placed on individual tracks. Due to a loss in tracking efficiency at low momentum, 
a lower transverse-momentum cut of 250 MeV/c was introduced for negatives. These 
conditions affect the A acceptance as shown in Fig. 2, middle. Additionally, a lower pt 
cut of 500 MeV/c for positives and an upper pt cut of 2.0 GeV/c for positives and 0.6 
GeV/c for negatives are introduced since they cost no signal and reduce background. 
Also, only pairs with an opening angle greater 20 mrad are accepted. For a clean A 
acceptance, we reject all A hyperons with pt < 0.9 GeV/c and pt > 2.5 GeV/c. There 
is also a fiducial cut in A rapidity: 2.0 < y < 2.4. Fig. 3 shows an Armenteros plot 
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Figure 4. Left: Invariant mass spectra for different extraction methods. On 
top are plotted the signal plus background (dashed line) and the normalized 
combinatorial background (solid line). The bottom shows the signal after 
background subtraction. 

Right: Number of reconstructed A as a function of the decay length. 



of the data and calculated lines for A and Kg. In order to isolate clean A particles, 
additional cuts were made on selecting 10 - 120 MeV/c and on the Podolanski- 
Armenteros variable a selecting values below 0.65, i.e. only A's where the proton is 
emitted backwards in the A rest frame. 

3. Consistency tests: 

In Fig. 4, left, the invariant mass spectrum is plotted for three different methods 
analyzing the same events; (1) TPC in a stand-alone mode, (2) TPC plus SDD veto 
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Figure 5. Left: Transverse momentum spectra for A shown separately for two 
different analysis methods and for the three centrality classes in the rapidity region 
2 < yA < 2.4. All spectra are acceptance corrected but not normalized to the 
number of events. 

Right: Midrapidity density dNyv/dy^jj as a function of the centrality of the 
collision. The plot shows data for the pt range of our measurement as well a an 
extrapolation to all pt (see text). 



and (3) TPC plus SDD veto plus vertex cut. One can notice a loss in signal but a 
large improvement in the signal-to-background ratio as successive cuts are introduced. 
The A signal, obtained after subtraction of combinatorial background, is fitted with 
a Gaussian. The fit parameters for all three methods are very similar, the average of 
mA — 1.117 ± 0.001 GeV/c'^ is close to the accepted value [5] and the resolution is tiA 
= 11.2 ± 0.4 MeV/c^ 

As an additional consistency test the acceptance corrected number of 
reconstructed A's is plotted as a function of the decay length in the A rest frame, 
as shown in the right part of Fig. 4. The superimposed solid line represents an 
exponential fit with a proper decay length ctq = (8.1 ± 0.5) cm, in good agreement 
with the accepted [5] value of 7.89 cm. 

4. Centrality classes 

The event sample was divided into three centrality classes according to the number of 
SDD1/SDD2 tracks as indicated in Tab. 1. Also indicated there are the corresponding 
fractions of the geometric cross section and the number of participants N^art evaluated 
using the UrQMD model [7] . 

5. Results 

The experimental data were corrected for the geometrical acceptance and 
reconstruction efficiency. In Fig. 5, left, the corrected but un-normalized A 
spectra are presented in different multiplicity bins for two different analysis 



Lambda production in 40 A GeV/c Pb-Au collisions 



6 



mult, class 


sddl 


sdd2 


sdd3 


SDD-tracks 


100 - 200 


200 - 300 


300 - 400 


events (10*^) 


0.8 


1.3 


0.7 


0-/CTgeo (%) 


15 - 36.7 


4.8 - 15 


< 4.8 


mean ATpart 


237 


328 


384 



Table 1. Definition of the three centrality classes 



methods. The invariant multiplicity is fitted with an exponential distribution 
dN/ptdpt oc exp{—mt/T\), also shown in Fig. 5. 
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Figure 6. Left: Dependence of the inverse slope parameter Ta on the centrality 
of the collision. In addition to the present data are shown data from AGS and full 
energy SPS as well as a calculation in the UrQMD model (references see text). 
Right; Mass dependence of inverse slope for the most central collisions. The 
CERES results for 40 AGcV/c (present data and [4]) are shown together with 
the systematics for full SPS energy [8] 

The dependence of the inverse slope parameter Ta on the number of participants 
Npart is presented in Fig. 6, left. The measured values of Ta exhibit a slight but 
not very significant increase as a function of centrality. As shown the slopes are 
close to those observed at top AGS and SPS energies. The measured values are 
not well reproduced by the UrQMD model (line in Fig. 6), as also observed for 
the slopes of the proton spectra [4]. The mass dependence of the inverse slope is 
presented in Fig. 6, right. Within errors, our present data agree with the systematic 
behavior of measurements at top SPS energies [8]. The large differences between the 
slopes (Th- = 176 ± 5 MeV, T(+)_(_) = 278 ± 12 MeV and Ta = 273 ± 20 MeV) 
indicate the presence of a strong radial flow also at 40 A GcV beam energy [4] . 

The midrapidity density dNA/dymid has been obtained by extrapolating the 
transverse momentum spectra to pt = using the functional form given above. On 
the right side of Fig. 5 the rapidity densities in the measured pt interval as well as 
extrapolated to pt = are presented as a function of centrality. A continuous rise 
with centrality is observed. In our most central multiplicity class we determine a total 
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midrapidity density dNA/dymid = 11.8 ± 2. A somewhat larger value of about 15 is 
reported by the NA49 collaboration [6]. Putting our present data into context with 
measurements at lower and higher beam energy, a continously rising behavior is seen 
as a function of -^s (see Fig. 7, left). 
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Figure 7. Left: dNA/dymitj as a function of the beam energy y/s (data for AGS 
and full energy SPS energies in [9, 10, 11]). 

Right: A/A ratio as a function of the beam energy y^s (other data from 
[12, 13, 11, 14]). 

We have also extracted the signal for the A. But due to the low statistics it is 
not possible to show transverse momentum spectra for A. Since all corrections are 
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Figure 8. Left: Comparison of data [17, 13] and particle ratios computed utilizing 
a thermal model [15]. 

Right: Phase diagram of nuclear matter [18] with the present data point for 40 
AGeV (square symbol). 
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symmetric for particle and antiparticle, they cancel in the ratio and we obtain Ra 
= A/A = 0.024 ± 0.010. The dominant error is statistical. Our result is in good 
agreement with the value recently reported by NA57 [13], albeit with a much larger 
statistical error for our present result. The A/A ratio as a function of -^/s (Fig. 7, 
right) shows the expected steep rise with increasing beam energies; for a net-baryon- 
free region the antiparticle/particle ratio is expected to be one and at RHIC energies 
this values is approached. 

Information about the chemical freozo-out parameters of the expanding system, 
the temperature T and the baryonic chemical potential /Ub , can be obtained from the 
ratios of particle yields. Utilizing the thermal model of [15] we have calculated several 
particle abundancies corresponding to the phenomenological condition [16] of fixed 
average energy per hadron of 1 GeV between top AGS and SPS energies, characterized 
by a parameterization in /Ub/T. The comparison of computed and measured particle 
ratios is presented in Fig. 8, left. Taking the different measured particle ratios into 
account we could extract the mean ratio /xb/T = 2.8±0.3. Following the Cleymans- 
Redlich freeze-out curve [16], this corresponds to values of /zb = 410 ± 30 MeV and T 
= 143 ± 5 MeV. This result is shown in the phase diagram Fig. 8, right. The chemical 
freeze-out point at 40 A GeV is near the phase boundary and falls into the systematics 
between top AGS and SPS energies. 
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